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ABSTRACT 

A recent analysis of a Lunar Laser Ranging (LLR) data record spanning 38.7 yr re- 
vealed an anomalous increase of the eccentricity e of the lunar orbit amounting to 
e m eas = (9 ± 3) x 10 -12 yr -1 . The present-day models of the dissipative phenomena 
occurring in the interiors of both the Earth and the Moon are not able to explain it. 
In this paper, we examine several dynamical effects, not modeled in the data analysis, 
in the framework of long-range modified models of gravity and of the standard New- 
tonian/Einsteinian paradigm. It turns out that none of them can accommodate e mea s- 
Many of them do not even induce long-term changes in e; other models do, instead, 
yield such an effect, but the resulting magnitudes arc in disagreement with e mcas - In 
particular, the general relativistic gravitomagnctic acceleration of the Moon due to 
the Earth's angular momentum has the right order of magnitude, but the resulting 
Lense-Thirring secular effect for the eccentricity vanishes. A potentially viable Newto- 
nian candidate would be a trans-Plutonian massive object (Planet X/Nemesis/Tyche) 
since it, actually, would affect e with a non-vanishing long-term variation. On the 
other hand, the values for the physical and orbital parameters of such a hypothetical 
body required to obtain at least the right order of magnitude for e arc completely 
unrealistic: suffices it to say that an Earth-sized planet would be at 30 au, while a 
jovian mass would be at 200 au. Thus, the issue of finding a satisfactorily explanation 
for the anomalous behavior of the Moon's eccentricity remains open. 



Key words: gravitation-Celestial mechanics-ephemerides-Moon-planets and satel- 
lites: general 



1 INTRODUCTION 



lAnderson fc Nietd i|2010l ), in a review of some astrometric 
anomalies recently detected in the solar system by several 
independent groups, mentioned also an anomalous secular 
increase of the eccentricitjQ e of the orbit of the Moon 



= (9 ±3) X 1CT 12 yr" 



(1) 



based on an analysis of a long LLR data record spanning 
38.7 yr (16 March 1970- 22 November 2008) performed by 
IWilliams fc Boggsl (|2009h with the suite of ac curate dynam- 
ical f o rce models of the DE 421 ephemerides (jFolkner et al.l 
l2008l : IWilliams et al.l l2008l ) including all relevant Newto- 
nian and Einsteinian effects. Notice that eq. (ft) is sta- 
tistically significant at a 3(7— level. The first presentation 
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1 It is a dimensionless numerical parameter for which < e < 1 
holds. It determines the shape of the Keplerian ellipse: e = 
corresponds to a circle, while values close to unity yield highly 
elongated orbits. 



of such an effect appeared in IWilliams et al.l |200ll ). in 
which an extensive discussion of the state-of-the-art in 
modeling the tidal dissi pation in both the Earth and the 
Moon was given. Lat e r. I Williams fc Dickevl (|2003l ). relying 
upon IWilliams et~ai] <|200ll ). yielded an anomalous eccen- 
tricity rate as larg e as e me as = (1.6 ± 0.5) x 10~ n yr -1 . 
lAnderson fc Nietoi |2010l) commented that eq. ([lj is not 
compatible with present, standard knowledge of dissipative 
processes in the interiors of both the Earth and Moon, which 
were, actually, modeled by IWilliams fc Boggsl (2009). Lhc 
relevant physical and osculating orbital parameters of the 
Earth and the Moon are reported in Table Q] 

In this paper we look for a possible candidate for ex- 
plaining such an anomaly in terms of both Newtonian and 
non-Newtonian gravitational dynamical effects, general rel- 
ativistic or not. 

To this aim, let us make the following, preliminary re- 
marks. Naive, dimensional evaluations of the effect caused 
on e by an additional anomalous acceleration A can be made 
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Table 1. Relevant physical and osculating orbital parameters of the Earth-Moon system, a is the scmimajor axis, e is the eccentricity. 
The inclination I is referred to the mean ecliptic at J2000.0. fl is the longitude of the ascending node and is referred to the mean equinox 
and ecliptic at J2000.0. Ul is the argument of pericenter. G is the Newtonian gravitational constant. The masses of the Earth and the 
Moon are M and m, respectively. The orbital parameters of the Moon were retrieved from the WEB interface HORIZONS (Author: J. 
Giorgini. Site Manager: D. K. Yeomans. Webmaster: A. B. Chamberlin), by JPL, NASA, at the epoch J2000.0. 



a (m) e / (deg) 


Q (deg) 


u (deg) 


GM (m 3 s" 2 ) 


m/M 


3.81219 x 10 s 0.0647 5.24 


123.98 


-51.86 


3.98600 x 10 14 


0.012 



by noticing that 



with 



A_ 

na ' 



na = 1.0 x 10 :! m s _1 = 3.2 x 10 10 m yr _1 



(2) 



(3) 



for the geocentric orbit of the Moon. In it, a is the orbital 
semimajor axis, while n = \J fi/a 3 is the Keplerian mean 
motion in which fi = GA/(1 + m/M) is the gravitational pa- 
rameter of the Earth-Moon system: G is the Newtonian con- 
stant of gravitation. It turns out that an extra-acceleration 
as large as 



A i 



3 x ICT 16 m s~ 2 = 0.3 m yr~ 2 



(4) 



would satisfy eq. JTJ). In fact, a mere order-of-magnitudc 
analysis based on eq. ((2} would be insufficient to draw mean- 
ingful conclusions: finding simply that this or that dynamical 
effect induces an extra-acceleration of the right order of mag- 
nitude may be highly misleading. Indeed, exact calculations 
of the secular variation of e caused by such putative promis- 
ing candidate extra-accelerations A must be performed with 
standard perturbative techniques in order to check if they, 
actually, cause an averaged non-zero change of the eccen- 
tricity. Moreover, also in such potentially favorable cases 
caution is still in order. Indeed, it may well happen, in prin- 
ciple, that the resulting analytical expression for (e) retains 
multiplicative factor^] l/e fe , k = 1, 2, 3, ... or e k , k = 1, 2, 3... 
which would notably alter the size of the found non-zero sec- 
ular change of the eccentricity with respect to the expected 
values according to eq. 

The plan of the paper is as follows. In Section [5] we deal 
with several long-range models of modified gravity. Section 
[jj] analyzes some dynamical effects in terms of the standard 
Newtonian/Einsteinian laws of gravitation. The conclusions 
are in Section [4] 



2 EXOTIC MODELS OF MODIFIED GRAVITY 
2.1 A Rindler-type acceleration 

As a practical example of the aforementioned caveat, let us 
consider the effective model for gravity of a cen tral object 
of ma ss M at large scales recently constructed bv lGrumillerl 
(|2010T ) . Among other things, it predicts the existence of a 
constant and uniform acceleration 



A = A 



Rinf 



(5) 



radially directed towards M. As shown in lloriol (|2010al ). 
the Earth-Moon range residuals 8p over At — 20 yr yield 
the following constrain for a terrestrial Rindler-type extra- 
acceleration 



Aitin < 5 x 10 -16 m s~ 2 = 0.5 m yr' 2 , 



(6) 

which is in good agreement with eq. Q. 

The problem is that, actually, a radial and constant 
acceleration like that of eq. ((5} does not induce any secu- 
lar variation of the eccentricity. Ind eed, from the stand ard 
Gaus£| perturbation equation for e (|Bertotti et al.|[2003l ) 



de 
lit 



Ar sin f + A T 



cos / + - (l - - 
e V a 



(7) 

in which / is the true anomahjj, and An, At are the radial 
and transverse components of the perturbing acceleration 
A, it turns out (jlorioll2010al ) 



Ae = 



Amn (l - e 2 ) (cosE - cos So) 



where E is the eccentric anomahjfl so that 



Ae| 2 " = 0. 



(8) 
(9) 



2.2 A Yukawa-type long-range modification of 
gravity 



It is well known that a variety of theoretical parad igms 
(|Adelberger et al.l 120031 : iBertolami fc Paramos! I2OO0D al- 
low for Yukawa-like deviations from the usual Newtonia n 
inverse-square law of gravitation (jBurgess fc Ck>utierlll988h . 
The Yukawa-type correction to the Newtonian gravitational 
potential ?7n = — where /j, = GM is the gravitational 
parameter of the central body which acts as source of the 
supposedly modified gravitational field, is 



U Y = 



afi 



' exp * ~ X 



(10) 



Here e denotes the eccentricity: it is not the Napier number. 



where is the gravitational parameter evaluated at dis- 
tances r much larger than the scale length X. 

In order to compute the long-term effects of eq. (| 10|) on 
the eccentricity of a test particle i t is convenient to ad opt 
the Lagrange perturbative scheme |Bertotti et al.|[2003l ). In 

3 It is just the case to remind that the Gauss perturbative equa- 
tions are valid for any kind of perturbing acceleration A, whatever 
its physical origin may be. 

4 It is an angle counted from the pericenter, i.e. the point of clos- 
est approach to the central body, which instantaneously reckons 
the position of the test particle along its Keplerian ellipse. 

5 Basically, E can be regarded as a paramctrization of the polar 
angle in the orbital plane. 



such a framework, the equa tion for the long-term variation 
of e is (|Bertotti et al.ll2003h 



de\ 
dt/ 



dTZ 



VI - e 2 Ouj 



dM ) 



(11) 



where to is the argument of pericenteiff], M = n(t — t p ) = 
E — e sin _E is the mean anomaly of the test particl^H, and 1Z 
denotes the average of the perturbing potential over one or- 
bital revolution. In the case of a Yukawa-type perturbatiorjf], 
eq. (fTU)l yields 



{Uy) = 



afioo exp I 



■f) r f ae 



(12) 



where Io(x) is the modified Bessel function of the first kind 
Ik(x) for k — 0. An inspection of eq. and eq. (| 12f) 

immediately tells us that there is no secular variation of e 
caused by an anomalous Yukawa-type perturbation which, 
thus, cannot explain eq. |T}. 



2.3 Other long-range exotic models of gravity 

The previous analysis has the merit of elucidating cer- 
tain general features pertaining to a vast category of 
long-range modified models of gravity. Indeed, eq. (|lip 
tells us that a long-term change of e occurs only if 
the averaged extra-potential considered explicitly de- 
pends on uj and on time through A4 or, equivalently, 
E. Actually, the anomalous potentials arising in the 
majority of long-range modified models of gravity are 



time- in dependent and sp h erical ly symmetric (IDvali et al 



2000; Capozzicllo et al 



2005 



2003 ; IDvali et al. 12003 



Gruzinov 



200ll; ICapozziello fc Lambiasd 



Kerr et al 



2003; 



... lAllemandi et al. 
120051 ; iJaekel fc Revnaudl l2005alE 



Apostolopoulos fc Te tradis 
ICapozziello et al 



2006a 




2006a 





IPC 



MoffatT 
Sanders ! 




Navarro fc van Acoleve nl 120051; iRevnaud fc Jaekell 12003: 



Brownstein fc Moffat 



Jaekel fc Revnaud 



Adkins et al. II2007I; iBertolami 



Capozzicllo 



Bertolami . . 

Ruggierd |2010| ; ISotiriou 



Francaviglia 2008 
Santos! 120091; Ide Fel 



Noiiri fc Odintsovl [2007 
TsuiikawaT I2OIC 



Navarro fc van Acoleven 
Adkins fc McDonnelll " 2007 
et al.ll2007l; ICapozziellol 2007 



Faraonil |2010| ; I Fabrina et al.l 



201 if ). Anomalous accelerations A exhibiting a dependence 
on the test particle's velocity v were also pr op osed in 
differe nt frameworks djaekel fc Revnaudl l2005al lbl; iHofaval 
l2009al lbl; iKehagias fc Sfetsoj l2009h . Since they have to be 



evaluated onto the unperturbed Keplerian ellipse, for which 



6 It is an angle in the orbital plane reckoning the position of the 
point of closest approach with respect to the line of the nodes 
which is the intersection of the orbital plane with the reference 
{x, y} plane. 

7 t p is the time of passage at pericenter. 

8 Several investigations of Yukawa-type effects on the lunar data, 
yielding more and more tight constr aints on its parameters, ar e 
present in the literature: see, e.g., iMiiller fe Biskupekl j2007Tl : 
liviiiller et al] J2007I . |2008| . l2009t) . 
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the following relations hold (|Murrav fc Correidliblch 



dt 



I'R 



Vt 



f: ' I i cos E) , 

nae sin E 



(13) 



1 — e cos E ' 



1 — e cos E ' 



where vr and vt are the unperturbed, Keplerian radial and 
transverse components of v, it was straightforward to infer 
from eq. (|7| that no long-term variations of the eccentricity 
arose at all (|loridl2007l ; llorio fc Ruggierdl2010l ). 

An example of time-dependent anomalous potentials 
occurs if either a secular change of the N ewtonian gravi- 
tational constant (|Milnelll935l ; |Diradll937h or of the mass 
of the central body is postulated, so that a percent time vari- 
ation of the gravitational parameter can be considered. 
In such a case, it was recently shown with the Gauss per- 
turbative scheme that the ecc entricity experiences a secular 
change given by (|loridl2010bh 

(e) = (l + e)(fy. (14) 

As remarked in llorid (|2010bh , eq. Q and eq. QD would 
imply an increase 

A 

If attributed to a change in G, eq. (|15|) would be one 
order of magnitude larger than the present-day bounds 
on GjG obtained f rorrC] LLR (IMiiller fc Biskupekl |2007l : 
I Williams et~al]|200^ . Moreover, iPitieval (|2010l T recently ob- 
tained a secular decrease of G as large as 



+8.5 x 10" 12 yr" 1 . 



(15) 



— = (-5.9 ±4.4) x 10 
G 



-14 -1 



(16) 



from planetary data analyses: if applied to eq. (|14|) . it is 
clearly insufficient to explain the empirical result of eq. ([l]). 
Putting aside a variation of G, the gravitational parameter 
of the Earth may experience a time variation because of 



(IBlinnikov & Khlopov 


19831; iKhlooov et al. 199ll; Khlopov 


1 19991; iFootJ 


20041; lAdler 


l2008t). iKhriplovich & Shepelvanskv 


(2009]) and 


Xu & Sicgcl (2008) assume for the Earth 



M 
M 



+10 yr , 



(17) 



which is far s maller than eq. (|15p . as noticed bv llorid (|2010d ). 
lAdlerl { 20081 ) yields an even smaller figure for M/M. 



3 STANDARD NEWTONIAN AND 

EINSTEINIAN DYNAMICAL EFFECTS 

In this Section we look at possible dynamical causes for eq. 
{TJ in terms of standard Newtonian and general relativistic 



9 According to iDirad lll937l i. G should decrease with the age of 
the Universe. 

10 Because of the secular tidal effects, the LLR-based determina- 
tions of G depend more strongly on the solar perturbations, and 
the (i/fi values should be interpreted as being sensitive to changes 
in the Sun's gravitational parameter GM. 
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gravitational effects which were not modeled in processing 
the LLR data. 



3.1 The general relativistic Lense-Thirring field 
and other stationary spin-dependent effects 

ft is interestin g to notice that the mag nitude of the gen- 
eral relativistic Lense fc Thirrinel (|l918T ) acceleration expe- 



rienced by the Moon because of the Earth's angular mom en- 
tum S = 5.86 x 10 33 kg m 2 s _1 (|McCarthv fc Petitll2004T ) is 
just 



2vGS 
c 2 a 3 



1.6 x 10~ 16 m s" 2 



0.16 m yr 



(18) 



i.e. close to eq. Q. On the other hand, it is well known that 
the Lense-Thirring effect does not cause long-term varia- 
tions of the eccentricity, fndeed, the integrated shift of e 
from an initial epoch c orrespondin g to fo to a generic time 
corresponding to / is (Soffel 1989) 



Ae = 



2GS cos I (cos / — cos fo) 
c 2 na 3 y/l — e 2 



(19) 



From eq. (|19[1 it straightforwardly follows that after one or- 
bital revolution, i.e. for / — > fo + 2n, the gravitomagnetic 
shift of e vanishes. In fact, eq. (|19p holds only for a specific 
orientation of S, which is assumed to be directed along the 
reference z axis; incidentally, let us remark that, in this case, 
the angle / in eq. (|19p is to be intended as the inclination 
of the Moo n's orbit with respect to the Earth's equatoj"] 
Actually, in lloriol |2010df ) it was shown that e does not sec- 
ularly change also for a generic orientation of S sine 



2Gn 



c 2 a(l 



[S z cos / + sin I (Sx sin Q — S y cos Q)] . 



(20) 

Thus, standard general relativistic gravitomagnetism cannot 

be t he cause of eq. flT] ). 

Iforio fc Ruggierol (|2009T ) explicitly worked out the grav- 
itomagnetic orbital effects induced on the trajectory of a test 
particle by the the weak-field approximation of the Kerr- 
de Sitter metric. No long-term variations for e occur. Also 
the general relativistic spin-spin effects a la Stern-G erlach 
do not cause long-term variations in the eccentricity (|loriol 
l2010dl ). 



3.2 General relativistic gravitomagnetic 
time-varying effects 

By using th e Gauss perturbative equations, 

iRuggiero fc Ioricl (|2010T ) analytically worked out the 
long-term variations of all the Keplerian orbital elements 
caused by general relativistic gravitomagnetic time-varying 



11 It approximately varies between 18 deg and 29 deg 
dSeidelmandll992l;lwilliams fc Dickey||2003h . 

12 Here /, fl can be thought as referring to the mean ecliptic 
at J2000.0. Generally speaking, the longitude of the ascending 
node Q is an angle in the reference {x,y} plane determining the 
position of the line of the nodes with respect to the reference x 
direction. 



effects. For the eccentricity, IRuggiero fc forid (|2010T ) found 
a non- vanishing secular change given by 



(e) 



GSi (2 + e)cos/ 



(21) 



in which 5*i denotes a linear change of the magnitude of the 

angular momentum of the ce ntral rotating bod y. 

in the case of the Earth, IRuggiero fc Ioriol (|2010T ) quote 



Si 



-5.6 x 10 16 kg m 2 s 2 



(22) 



due to the secular decrease of the Earth's diurnal rotation 
period (|f3rosche fc Schuhlll998T l P/P = -3 X 10~ 10 yr" 1 . 
Thus, eq. (|21[) and eq. (|22|l yield for the Moon's eccentricity 

(e) = -2 x 10~ 23 yr~\ (23) 

which is totally negligible with respect to eq. Q. 



3.3 The first and second post-Newtonian static 
components of the gravitational field 

Also the first post-Newtonian, Schwarzschild-type, spheri- 
cally symmetric static component of the gravitational field 
which was, in fact, fully modeled by IWilliams fc Boggj 



does not induce long-term variations of e ((Soffel 
The same holds also for the spherically sym- 
metric second pos t-Newtonian terms of order Q(c 



(|Damour fc Schaferl Il988l ; ISchafer fc Wexlll993l; IWexlll995l ) , 
which were not modeled by I Williams fc Boggsi (|2009h . fn- 
deed, let us recall that the components of the spacetime 
metric tensor <? M „, n, v = 0, 1, 2, 3, are , up to the second 
post-Newtonian order, (|Nordtvedtll 1996T I 



floo 



+ 21 



9ij 



!(f) 2 



+ 



- -^[l + 2f + | (f )+...], = 1,2,3, 

(24) 

where 3JI = p/c 2 . Notice that eq. (|24p are written in the 
standard isotropic gauge, suitable for a direct comparison 
with the observations. Incidentally, let us remark that the 
second post-Newtonian acceleration for the Moon is just 



2 2 

fj, n 



4x 10" 25 m s" 2 



4x 10 -10 m yr" 2 . (25) 



3.4 The general relativistic effects for an oblate 
body 

ISoffel et al.l l| 19881 ). by using the Gauss perturbative scheme 
and the usual Keplerian orbital elements, analytically 
worked out the first-order post-Newtonian orbital effects in 
the field of an oblate body with adimensional quadrupole 
mass moment J2 and equatorial radius R. 

It turns out that the eccentricity undergoes a non- 
vanishing har monic long-term v ariation which, in general 
relativity, if0 (|Soffel et al.lll988T l 



13 Here ui refers to the Earth' s equator, so that its period 
amounts to 8.85 yr jRoncolill2005l 'l. 
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In view of the fact that, for the Earth, it is J2 = 
1.0826 3 x 10~ 3 dMcCarthv fc Petitl l 2004h and R = 6.378 x 
10 3 m (jMcCarthv fc Petidl2004l ). it turns out that the first- 
order general relativistic Ji,c~ 2 effect is not capable to ex- 
plain eq. |T} since it is 



(e) <4 x 10 



-19 -1 

yr 



(27) 



as a limiting val u e for the periodic perturbation of eq. (|26|) . 

ISoffel et al.l (|l988l ) pointed out that the second-order 
mixed perturbations due to the Newtonian quadrupole 
field and the general relativistic Schwarzschild accelera- 
tion are of the same order of magnitude of the first-order 
ones: their orbita l effect s were analytically worked out by 
iHeimberger et alj (|l99Ch with the technique of the canoni- 
cal Lie transformations applied to the Delaunay variables. 
Given their negligible magnitude, we do not further deal 
with them. 



3.5 A massive ring of minor bodies 

A Newtonian effect which was not modeled is the action 
of the Trans -Neptunian Objects (TNQs) o f the Edgeworth- 
Kuiper belt (|Edgeworthll 19431 ; iKuiperil 195 lb ■ It can be taken 
into account by means o f a massive cir cular ring having mass 
m ri n g < 5.26xl0~ 8 Mq jPitieval201Ch and r adius R ling = 43 
au jPitieval [20ich ■ Following iFienga et all (|200Sft , it causes 
a perturbing radial acceleration 



A r ; ne — 



2 



2rR. 



[b ( 3 >{a) -a6 ( 3°'(a)] r, a 



Rring 



(28) 



The Laplace coefficients are defined as jMurrav fc Dermottl 
1 19991 ) 



= - 

n 



cos(jV') 



-dip, 



(29) 



(1 - 2acos?/' + a 2 ) s 

where s is a half-integer. Since for the Moon a ~ 3 x 10 -10 , 
eq. (f28|) becomes 



with 



-*4-rin 



10" 



ar, 



10" 



m yr 



(30) 



(31) 



which is far smaller than eq. (j4]). 

Actually, the previous results holds, strictly speaking, 
in a heliocentric frame since the distribution of the TNOs is 
assumed to be circular with respect to the Sun. Thus, it may 
be argued that a rigorous geocentric calculation should take 
into account for the non-exact circularity of the TNOs belt 
with respect to the Earth. Anyway, in view of the distances 
involved, such departures from azimuthal symmetry would 
plausibly display as small corrections to the main term of 
eq. (|28|l . Given the negligible orders of magnitude involved 
by eq. (|3ip . we feel it is unnecessary to perform such further 
calculations. 

The dynamical act ion of the belt of the minor asteroids 
l|Krasinskv et al.ll2002T ) was, actually, modeled, so that we 
do not consider it here. 



3.6 A distant massive object: Planet 
X /Nemesis /Tyche 

A promising candidate for explaining the anomalous 
increase of the lunar eccentricity may be, at least in 
principle, a trans-Plutonian massive body of plane- 
tary size located in the remote perip heries of the solar 
syste m: Planet X/Nemes is /Tych e dLvkawka fc Mukail 



200 



2011 



.Melott fc Bambachl |2010| ; iFernandej 
Matese fc WhitmireT I 201 ll ). Indeed, as we will see, the 
perturbation induced by it would actually cause a non- 
vanishing long-term variation of e. Moreover, since it 
depends on the spatial position of X in the sky and on its 
tidal parameter 



ICx = 



Gmx 



(32) 



where mx and dx are the mass and the distance of X, respec- 
tively, it may happen that a suitable combination of them 
is able to reproduce the empirical result of eq. |T}. 

Let us recall that, in general, the perturbing poten- 
tial felt by a test particle orbiting a central body due to 
a very distant, poi ntlike mass can b e cast into the following 
quadrupolar form (|Hogg et al.l ll99iT) 



Ux = 



3 [r-l 



(33) 



where I = {l x ,l y ,lz} is a unit vector directed towards X 
determining its position in the sky; its components are not 
independent since the constraint 



~i~ ly ~H ^2 



(34) 



holds. By introducing the ecliptic latitude /3x and longitude 
Ax in a geocentric ecliptic frame, it is possible to write 



l x = cos /3x cos Ax, 



cos /3x sin A> 



sin/3 x . 



(35) 



In eq. (f33|) r — {x, y, z} is the geocentric position vector 
of the pertu r bed p article, which, in the present case, is the 
Moon. Ioriol (|201ll ) has recently shown that the average of 
eq. (|33[1 over one orbital revolution of the particle is 



(t/x) 



K,xi' 
32 



U 



(36) 



with U ^e, /, £1, u; tj given by eq. (|37[). Note that eq. (|36[) 
and eq. (|37[) are exact: no approximations in e were used. In 
the integration I was kept fixed over one orbital revolution 
of the Moon, as it is reasonable given the assumed large 
distance of X with respect to it. 

The Lagran ge planetar y equation of eq. ()11[) straight- 
forwardly yields l)loridl201lh 



15K,xeVT~e5 
16n 



£ I,Sl,u):l 



(38) 



with £ (j, fl, oj; tj given by eq. (J39J) . 

Actually, the expectations concerning X are doomed to 
fade away. Indeed, apart from the modulation introduced by 
the presence of the time- varying /, uj and Q, in eq. (|39|l , the 
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U = — (2 + 3e 2 ) (—8 + 91% + 9l y + 6Z 2 ) — 120e 2 sin2aj (l x cos Q + l y sinQ) [l z sin 7+ 

+ cos/ (ly cosfi - l x sinf2)] - 15e 2 cos 2^ [3 (l x - l 2 y ) cos2^ + 2 (l x + Z 2 - 2Z 2 ) sin 2 I- 

— 4Z Z sin 27 (l y cos f2 — sin U) + 6ZzZ H sin2Q] - 6 (2 + 3e 2 ) [(Z 2 - Z 2 ) cos 2Q sin /+ 

+ 2Z Z sin 2/ (i H cosQ - Z^ sinfi) + 2^;,, sin 2 7 sin 2f2] - 3 cos 21 {(2 + 3e 2 ) (Z 2 + Z 2 - 2Z 2 ) 

+ 5e 2 cos2w [(Z 2 - Z 2 ) cos2Q + 2;^ sin2f2] } . 



(37) 



= — 8l z cos 2lj sin I (l x cos f2 + l y sin J7) + 4 cos J cos 2oj [—2l x l y cos 2f2+ 

+ (Z 2 - Z 2 ) sin2H] +sin2w [(Z 2 - Z 2 ) (3 + cos 21) cos 2Q + 2 (Z 2 + Z 2 - 2Z 2 ) sin 2 /- 

— 4Z Z sin 2/ (Zy cos £1 — l x sin f2) + 2l x l y (3 + cos 27) sin 2f2] . 



(39) 



values for the tidal parameter which would allow to obtain 
eq. CO are too large for all the conceivable positions {/3x, Ax} 
of X in the sky. This can easily be checked by keeping oj and 
Q. fixed at their J2000.0 values as a first approximation. 

Figure [T] depicts the X-induced variation of the lunar 
eccentricity, normalized to eq. CO. as a function of /3x and 
Ax for the scenarios bv lLvkawka fc Mukail J200S|) fm? ax = 
0.7 m m ,dT n = 101.3 au). and by iMatese fc Whitmird 
|201ll ) (mx" = 4 mj up ,dx = 30 kau). It can be noticed 
that the physical and orbital features of X postulated by 
such two recent theoretical models would induce long-term 
variations of the lunar eccentricity much smaller than eq. 
{TJ. Conversely, it turns out that a tidal parameter as large 
a.s 



£x = 4.46 x 10" 



(40) 



would yield the result of eq. CO. Actually, eq. (|40p is to- 
tally unacceptable since it corresponds to distances of X as 
absurdly small as dx = 30 au for a terre strial body, and 
dx — 200 au for a Jovian mass (|loridl201ll ). 

We must conclude that not even the hypothesis of 
Planet X is a viable one to explain the anomalous increase 
of the lunar eccentricity of eq. CO. 



4 SUMMARY AND CONCLUSIONS 

In this paper we dealt with the anomalous increase of the 
eccentricity e of the orbit of the Moon recently reported from 
an analysis of a multidecadal record of LLR data points. 

We looked for possible explanations in terms of unmod- 
eled dynamical features of motion within either the stan- 
dard Newtonian/Einsteinian paradigm or several long-range 
models of modified gravity. As a general rule, we, first, no- 
ticed that it would be misleading to simply find the right 
order of magnitude for the extra-acceleration due to this or 
that candidate effect. Indeed, it is mandatory to explicitly 
check if a potentially viable candidate does actually induce 
a non-vanishing averaged variation of the eccentricity. This 
holds, in principle, for the search of an explanation of any 
other possible anomalous effect. Quite generally, it turned 
out that any time-independent and spherically symmetric 
perturbation does not affect the eccentricity with long-term 
changes. 

Thus, most of the long-range modified models of gravity 



proposed in more or less recent times for other scopes are 
automatically ruled out. The present-day limits on the mag- 
nitude of a terrestrial Rindler-type perturbing acceleration 
are of the right order of magnitude, but it does not secularly 
affect e. As time-dependent candidates capable to cause sec- 
ular shifts of e, we considered the possible variation of the 
Earth's gravitational parameter fi both because of a tem- 
poral variation of the Newtonian constant of gravitation G 
and of its mass itself due to a steady mass accretion of non- 
annihilating Dark Matter. In both cases, the resulting time 
variations of e are too small by several orders of magnitude. 

Moving to standard general relativity, we found that 
the gravitomagnetic Lense-Thirring lunar acceleration due 
to the Earth's angular momentum, not modeled in the data 
analysis, has the right order of magnitude, but it, actually, 
does not induce secular variations of e. The same holds also 
for other general relativistic spin-dependent effects. Con- 
versely, e undergoes long-term changes caused by the general 
relativistic first-order effects due to the Earth's oblateness, 
but they are far too small. The second-order post-Newtonian 
part of the gravitational field does not affect the eccentricity. 

Within the Newtonian framework, we considered the 
action of an almost circular massive ring modeling the 
Edgeworth-Kuiper belt of Trans-Neptunian Objects, but it 
does not induce secular variations of e. In principle, a vi- 
able candidate would be a putative trans-Plutonian massive 
object (PlanetX/Nemesis/Tyche), recently revamped to ac- 
commodate certain features of the architecture of the Kuiper 
belt and of the distribution of the comets in the Oort cloud, 
since it would cause a non-vanishing long-term variation of 
the eccentricity. Actually, the values for its mass and dis- 
tance needed to explain the empirically determined increase 
of the lunar eccentricity would be highly unrealistic and in 
contrast with the most recent viable theoretical scenarios for 
the existence of such a body. For example, a terrestrial-sized 
body should be located at just 30 au, while an object with 
the mass of Jupiter should be at 200 au. 

Thus, in conclusion, the issue of finding a satisfactorily 
explanation of the observed orbital anomaly of the Moon 
still remains open. Our analysis should have effectively re- 
stricted the field of possible explanations, indirectly pointing 
towards either non-gravitational, mundane effects or some 
artifacts in the data processing. Further data analyses, hope- 
fully performed by independent teams, should help in shed- 
ding further light on such an astrometric anomaly. 
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Figure 1. Long-term variation of the lunar eccentricity, normalized to cq. (T), induced by a trans-Plutonian, pointlike object X as a 
function of its ecliptic latitude /3x and longitude Ax- The node H and the perigee u> of the Mo on were kept fixed t o the J2000.0 values 
quote d in Table[T] The scenarios for the perturbing body X are those bv lLvkawka fc M ukai ( 200<|) (left panel), and bv lMatese fc Whitmird 
1120111 1 (right panel). 
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